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Abstract 

A Kallen-Lehman approach to 3D Ising model is analyzed numerically 
both at low and high temperature. It is shown that, even assuming a min- 
imal duality breaking, one can fix three parameters of the model to get a 
very good agreement with the MonteCarlo results at high temperatures. 
With the same parameters the agreement is satisfactory both at low and 
near critical temperatures. How to improve the agreement with Monte- 
Carlo results by introducing a more general duality breaking is shortly 
discussed. 
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Introduction 

The three dimensional Ising model (henceforth 3DI) is one of the main open 
problems in field theory and statistical mechanics. A large number of interest- 
ing statistical systems near the transition point are described by 3DI and the 
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theoretical methods suitable to deal with such a problem manifest deep connec- 
tions in various areas of physics ranging from quantum information theory to 
string theory (for a review see e.g. [T]). Besides its intrinsic interest in statistical 
physics, since the formulation of the Svetitsky-Yaffe conjecture [5], it has been 
widely recognized its role in describing the deconfinement transition in QCD. 
For this reason, the 3D Ising model is worth to be further investigated. 

It has been recently proposed ^ to use in statistical mechanics the powerful 
tools of Regge theory which have been so fruitful in the study of the strong 
interactions leading to the formulation of the dual models [5] (two detailed 
reviews are [7] and [H])- It has been argued in [S] that such ideas may be useful 
in dealing with the 3D Ising model. It would be also interesting to try to 
test these methods in the mean field theory of spin glassQ as developed in [13] 
(recently it has been proved, that such a framework provides one with the exact 
solution in the mean field Ising spin glass [l4l [T5l fT6| ) . 

In the present paper we perform a numerical test of these methods in the 
case of the 3D Ising model. We will try to find the best "Regge parameters" 
appearing in the Kallen-Lehman form for the free energy: that is the parame- 
ters which provide one with the best possible agreement both at high and low 
temperatures MonteCarlo expansions as well as around the critical point. 

The paper is organized as follows: in the second section the Kallen-Lehman 
form for the free energy of the 3D Ising model is discussed. In the third sec- 
tion the issue of duality breaking is analyzed. In the fourth and fifth sections 
the comparison with MonteCarlo results at high and low temperatures are pre- 
sented. In the sixth section, the behavior at the critical point is described. In 
the seventh section some possible improvements are pointed out. Eventually, 
we draw some conclusions. 

^Detailed reviews on the subject are |10) |12) . 
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1 The Kallen- Lehman form for the free energy 



The Kallen-Lehman representation [9] gives rise to an ansatz for the free energy 
of the 3DI model of the following fornix 
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here, for sake of simplicity, wc consider C^i — 1/2 = ^3 (the 2D Ising values). In 
order to avoid a complex free energy (which could arise due to the two square 
roots). 
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At high temperatures it is convenient the variable t = tanh (3 while u = exp(— 2/3) 
is the standard variable at small temperatures. Thus our trial free energy will 
depend on 3 parameters A, C2 ^7 Co C to be fixed in comparison with 
MonteCarlo data. 

The natural value (suggested by the application of Regge theory) of the 
factor q in the trivial one dimensional term in Eq. ([3]) is g = 2 (and this is 
the value that we will consider from now on, see |9j) since the one dimensional 
term of exact solution of the 2D Ising model is precisely log 2 cosh 2/3. On the 
other hand, we have also obtained quite good numerical results for q = \: in 
particular, for the set of parameters 



C* 1 = 1.096, v*^ = 2.586, A* 1 = 0.127 



(7) 



^The "Regge parameters" appearing in tliese formulas are related with those appearing in 
[9] by the following identities C,-^ = u, C,^ = 1/2 = <^2i Co ~ 
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we get the following deviations at high and small temperatures 



'^HT{(,q=l,l'q=lAq=l) ~ "W ^ ~W (8) 
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(a definition of tht and ult will be provided in the next sections in Eqs. (|20p . 
((22|) . ((25| and (US])). These quite small deviations obtained with q = 1 (which 
indeed is less natural than the ansatz with q = 2) show that the present method 
is well suited to describe the thermodynamics of the 3D Ising model. In any 
case, as it will be shown in a moment, the results for q = 2 are quite better. 

Besides the first trivial term in Eq. ([3]), the "Regge free energy" depends 
on the temperature through an effective fce//(/3) in the same way as the (non 
trivial part of the) free energy of the 2D Ising model depends on /3: 

(k ^ V / exp(2/3) - exp(-2/3) V 

- Uosh2/3coth2/3j ' ^ ^ 

The duality symmetry of the 2D Ising model (discovered in |17j by Kramers 
and Wannicr before the exact solution of Onsager [TB]) is manifest in the above 
function of the temperature since fc(/3) has the same expression when rephrased 
in terms of the low temperatures variable u and the high temperatures variable 
t. Indeed, the 3D Ising model has not such a duality symmetry: in the next 
section it will be described how the simplest duality breaking can be achieved. 



2 Minimal duality breaking 

There are indeed many reasonable ways to break duality, how can one choose? 
Here a criterion of simplicity will be followed. An important ingredient is the 
following: all the Ising models on simple hypercubic lattices are invariant under 
the Marchesini-Shrock transformation ^19j 

TT 

/3 — > /3 + ^T^-^i n G Z. 

While this symmetry is obviously realized in the 2D case in which one has at own 
disposal the exact solution, such an exact symmetry puts a strong constraint 
on the form of the free energy. As discussed in 9J , the free energy suggested by 
the Kallen-Lehman representation is expect to depend on (a suitable) fce//(/3). 
A simplicity criterion (which, as it will be shown in a moment, is supported by 
the numerical data) suggests to consider the easiest possible modification of the 
k of the 2D case in Eq. (fTO| : 

4 da exp(2/3) - ^2 exp(-2/3) 
di (exp(2/3) + rfocxp(-2/3))^ 



keff[P) = —- , , , ^^,,2 ' (11) 
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the fact that the exponential only depend on 2/3 is the simplest possible way to 
fulfil the Marchesini-Shrock symmetry. The four parameters di {i = 0, ..,3) are 
not independent since there are two conditions that the function kef / has to ful- 
fil: firstly, the maximum of k^f / has to occur at the known critical temperature 
(/?*) ^ of the 3D Ising model: 

f3*\ df3keffm^=p' = (12) 
/?* = 0.22165. (13) 

The second condition is related to the fact that the expected transition has to 
occur [9] when keff{P) = 1 as it happens in the 2D case: 

keffin = 1. (14) 

Thus, in order to consider the simplest modification of the two-dimensional fc(/3) 
in Eq. PTI)) . we will take = 1 (the 2D Ising value) in such a way that only one 
parameter is left: di and d2 can be expressed in terms of do and j3* as follows 

^ _ exp(4/9*) - 3do _ ^^^^ 



di 



3-rfo exp(-4/3*)' 

(exp(2^*) + doexp(-2/3*))' 



(16) 



From the numerical point of view, this parametrization is very useful since auto- 
matically ensures that the maximum is equal to one avoiding possible problems 
related to negative numbers appearing inside the square root in Eq. ([2]). More- 
over the minimal modification criterion suggests that 

do = 1 (17) 

as in the 2D case. 

Indeed, one should expect a more general duality breaking: the 2D Ising 
model and its exact duality are closely related to N=4 SUSY Yang-Mills theory 
(for instance see [SO])- From the gauge theory side, one would thus expect that 
a minimal duality breaking should be related to a Yang-Mills theory with two 
supersymmetry which has an effective duality (for instance see |21jV The 3D 
Ising model is related to QCD without supersymmetries, so that it is natural to 
expect that a more general pattern of duality breaking (which will be discussed 
later on) should occur in the 3D Ising model. However, it is a truly remarkable 
feature of the present tools that already a minimal duality breaking works very 
well in comparisons with MonteCarlo results. It is worth to mention here that 
the so called " effective string approch" to the 3D Ising model (in which duality 
symmetry is unbroken) gives reasonable resuts in comparison with MonteCarlo 
data (see for instance f22]). 
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One can express fee// in terms of the high and low temperatures variables: 

^e//H = -T (18) 
"1 (1 + dow ) 



. 4(1 -d,)- 2,(1 + -i^j + d-i.)] 

lii [P(l+<i„)-2((l-<io) + (l+<io)'' 



so our ansatz for the 3D Ising free energy is given by Eq. ([T]) in which k^ff is 
in Eq. ([TT]) with the constants in Eqs. HH) and p7)) . 

3 High temperatures 

The idea is to find the best set of high temperatures parameters (Cq, C2J •^*) in 
Eqs. ((H), dSl) (with the constraints in Eqs. ([13, (US]) and (HH)) which repro- 
duces as close as possible the available Monte Carlo data (see [23]). A hypercubic 
lattice has been chosen in the parameters space (every point in the lattice rep- 
resenting a possible set of high temperature parameters) , then the free energy 
^ will be evaluated at every point in the lattice. The best choice of parame- 
ters will be the one minimizing the following deviation function which, to some 
extent, represents the deviation between the ansatz and the Monte Carlo data: 
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P^-P^-l = ^ , /350 = /3ma. = 0.03, /q - 2.4831 

iP) = 3cosh/3-|- (^3(tanh/3)^ + 22(tanh/3)Vl87.5(tanh/3)^ 
1980 (tanh/3)^° + 24044 (tanh/3)^^ + 319170 (tanh/3)^'' 



where we keep the terms up to the 15th of [23j since our algoritm is not sensitive 
to the higher order terms, F iJ^'^ is the high temperatures MonteCarlo free 
energy, (3^ can be assumed to be of order O.O^O and /q is a constant introduced 
for numerical convenienccQ. 

In order to compare the Regge coefficients with the high temperature coeffi- 
cients in [23], we need to change variable from (3 to t = tanh/?, so that we have 
to use the expression in Eq. (|19p whenever it appears kef / . 



After /3 0.05 one is not anymore in the high temperatures regime since the critical 
temperature is at /3* 0.22 which is only a factor of four larger. Indeed, f) 0.03 appears 
to be not enough smaller than /3* . Nevertheless, we will see that the agreement of the Regge 
free energy with MonteCarlo data is excellent up to /3 fs; 0.03. 

■'in order for the program to be able to select good points in the parameters space an 
arbitrary constant has to be added. Otherwise, the C + + program would select uncorrect 
functions: for instance, it would select two curves which intersect in a point near /3 = (but 
whose shapes are very different) instead of two parallel curves. 
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The set of points candidates to be the best parameters are: 

C = 1.9389, u* = 1.9205, = 0.0781 (21) 

{2'Ky 

In order to have a graphical idea of the closeness of the Regge and the 
MonteCarlo free energies at high temperature we have drawn in Fig. ([T]) both 
functions. The agreement appears to be exceUent: the deviation at high tem- 
perature is 

a^T(r,^*,A*)«yS^^^^^«10-«. (22) 

It is worth to note that in the range of parameters 

1 < C* < 2, 1.5 < ly* < 2.5, 0.05 < A* < 0.2 

we found many other good sets of parameters with only shghtly larger deviations 
with respect to the MonteCarlo free energy. However, as it will be discussed in 
a moment, it is very difficult to evaluate the errors on the single coefficients of 
the expansions. Thus, we simply considered the set in Eq. (PT|) but other sets 
of parameters in the above range could reproduce better the single coefficients. 

Figure 1: MonteCarlo F^jj^'^'-^ (j3) and Kallen-Lehman i^g^ \f3) free ener- 
gies at high temperatures versus p. 

3.0014 I 1 1 , 1 , 1 




There are some points in the parameters space in which the agreement at low 
temperatures is much better than the agreement at high temperatures (which is 
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in any case quite good). In the present analysis it is apparent that one cannot 
achieve an excellent agreement both at high and at low temperatures: in the 
best cases very small deviations (of the order of 10~^) on one side correspond 
on the other side to quite good but not excellent deviation of the order of 10""*. 
This is a manifestation of the already discussed fact that the 3D Ising model, 
being related to QCD without any supersymmetry, needs a more general duality 
breaking. 



3.1 Numerical Derivatives at high temperatures 

A consistency test would be provided comparing the derivatives of MonteCarlo 
free energy with the ones predicted by our model at /? = 0. Due to the com- 
plicated expression of the Regge free energy numerical methods are needed. To 
this aim we consider the finite difference approximation of the derivative which, 
at any order of derivation n, is given by the forward finite difference of order n 
divided by the increase of the independent variable (say 5) raised to the n-th 
power, i.e. 

A^[/](a;) being the n-th order finite difference. Behind this apparently easy 
procedure there are hidden dangers which already at the lowest orders deserve 
special attention (for an introductory discussion see [25 )• which is the best 
choice of 5 in order to approximate correctly the value of the n-th order deriva- 
tive (since neither too small nor too large values are correct)? One can only try 
to find the suitable ^(„) to reproduce the MonteCarlo coefficients at any order. 
In the case of the Regge free energy it's always possible to find (more than) one 
suitable numerical value for 5. One can verify that the (5(„) which reproduce the 
MonteCarlo results fulfil a sort of "natural" scaling (see [21]): namely, as the 
order of derivation increases, the suitable 5 reproducing the MonteCarlo result 
stays constant or decreases. Despite being only a first attempt to compute the 
derivatives of the Regge partition function in zero, nevertheless it is a quite en- 
couraging fact that at any order one can find at least one 5 reproducing exactly 
the MonteCarlo data. 

Another possible way to obtain the numerical derivatives for /3 = is the fol- 
lowing. One can fit the Regge free energy (in which now the parameters are fixed 
to be the best parameters) for small /3 with a polynomial with unknown coefh- 
cients (in which a priori one has to include all the-even and odd-coefficients) in 
such a way that the coefficients of the polynomial are the numerical derivatives 
at the origin. However, after the first four coefficients are fixed, the polynomial 
becomes very close to the Regge free energy. The higher order coefficients are 
fixed by the program to be the MonteCarlo results, but it is very difhcult to 
evaluate the errors since if one changes the higher order coefficients of a factor of 
two or four the deviation with respect to the optimal polynomial is smaller than 
10~^. Even with this method, it is extremely difficult to evaluate the deviations 
of the derivatives of the higher order derivatives of the Regge free energy from 
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the MonteCarlo results. In any case, we find an attractive feature of our model 
that even fixing the first few coefficients at high temperature one obtains a good 
agreement at the critical point and at low temperatures. 



4 Low temperatures 

Once the best set of high temperatures parameters {(* , v* , X*) have been found 
in Eq. ([2T|) one could hope to find a good agreement also at low temperatures 
(something which, a priori, is very far from obvious). The internal energy at 
low temperatures is (fee// as in Eq. (fT8|) and expressing Fm as in ([3|) in terms 
of u) 



\ {u) + 2h^2u^F^'^j^-'^Hu) . (23) 



KL 

This is the average energy for spin (the two expressions can differ by a non zero 
constant /i, see footnote 4) to be compared with {■^}^^^ (the polynomial in u 
representing the MonteCarlo average energy for spin for small u) found in [25] . 

It is worth to discuss the origin of the constant Ii in the above formula. The 
partition function of the 3D Ising model reads 

Z^D = ^exp(-/3(i/o)), 
Hq = aiffj + const 

(ij) 

where the constant in the Hamiltonian can always be chosen in such a way that 
the lowest energy state has zero energy. The ansatz ^ has not been deduced 
by solving in some approximated way the model. Rather, it has been found 
on the basis of the Regge theory and the Kallen-Lehman representation: one 
cannot pretend that it already corresponds to the normalization in which the 
Hamiltonian Hq has a zero energy ground state. At very low temperatures, the 
free energies corresponding to Hq and to Hq+ const differ by a term proportional 
to l/T 



log 



exp (— /3 (Hq + const)) 



-log 



^exp(-/3(i7o)) 



(3r^ l/T~logu. 



In terms of the internal energy, such a term corresponds to a constant: 

2u— logu = const, 
au 

From the numerical point of view it is extremely inconvenient to use the 
expression on the right hand side of Eq. ((23|) since the derivative du has a 
really big expression. A trivial but useful trick is the following: in [25j one has a 
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polynomial expression for the average energy for spin so that one can integrate it 
and obtain the Monte Carlo expression for the free energy at low temperatures: 



P\KL ^ 14 



this implies that 



"Mm) 



li* — 2/1 log u 



Therefore, the low temperatures test function reads 

50 

1=1 



Ui — Urain^ ^50 ~ ^max — 0.3, 



(24) 
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+ h log u 



where Ii = —1 and I2 = 0.1034. We keep the terms up to the Ibth of [25] since 
our algoritm is not sensitive to the higher order terms and Umax can be assumed 
to be of order 0.3. 

Actually, the value of u which corresponds to the critical temperature is 
Ucrit = crit) ~ O-^- The low tcmpcratures expansion provides one with 

reasonable results for u <C 0.6. Indeed, Umax ~ 0.3 is not at all much smaller 
than Ucrit- one may expect a reasonable value of Umax to be quite smaller tharj^l 
0.01. Nevertheless, we will see that the agreement is quite good up to u w 0.3. 

Eventually, the deviation at low temperatures between the Regge and the 
MonteCarlo free energies evaluated for the same optimal parameters in Eq. (|21[) , 
which have been found by asking the optimal agreement at high temperatures 
is 



^/Xlt{C,v*,\*) 



50 



10 



-5 



(26) 



and one can see that the agreement at low temperatures is two order of magni- 
tude better than in the case in which in Eq. ^ q = 1 (see Eq. ©). 



^Unfortunately, a numerical analysis in the region of m < 0.01 is quite difficult because 
of the fact that the logaritm is very large and dominates the too small numbers coming 
from the MonteCarlo polynomial. If the logartmic divergent terms would be absent, in the 
region < w < 0.01 very soon the MonteCarlo polynomial would give rise to very small 
numbers (below the precision of our PC). Thus, to analyze this range more refined numerical 
techniques are needed. In any case, the good agreement at the critical point provides one with 
some control on the deviation of the theoretical low temperatures coefficients with respect to 
the MonteCarlo ones. 
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Figure 2: MonteCarlo F^j,^^^ and Kallcn-Lchman F^^j^''"' •^'\u) free energies 
at low temperatures versus u. 



"LTPoLMonteCarlo" 
"LTJheor_KL" 



Remarkably enough, at low temperatures also the agreement is quite good 
as can be also seen in Fig. ([2]). As far as the numerical derivatives at w = 
of Regge free energy, the same considerations as in the high temperatures case 
hold. Namely, both a direct evaluation by means of a simple algoritm and to 
try to find the suitable approximating polynomial give results compatible with 
the MonteCarlo ones but it is not clear to the present authors how to estimate 
the errors on the single coefficients. It is nevertheless worth to stress here that 
the relative deviations both at high and at low temperatures are very small (see 
Eqs. ([22|) and (|26|)). Furthermore, as it will be now shown, also the critical 
exponent is in a very good agreement with the available data: one may hope 
that the good agreement at the critical point could prevent too large deviations 
of the the Kallen-Lehman derivatives from the MonteCarlo data. 



5 The Critical Point 

The last test is the critical point. Namely, once the parameters have been fixed 
as in Eq. (|2ip one may hope to verify that the critical point is correctly predicted 
too. To do this one can fit near the critical point the non-analytic part of the 
free energy in Eq. with the optimal parameters with a function of the form 

^That is, one has to exclude the term log 2 cosh /3 which does not contribute to the critical 
behavior. 
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FcRIT 

FcRiT~c,\(3-f3*f-'' + C2 (27) 

(where ci, C2 and a are constant) and find the optimal values for ci, C2 and 
a so that a will be our estimate for the critical exponent. This form of the 
free energy's critical part is expected both from Conformal Field Theory and 
experiments. The results of the fit done with the PC program Mathematica 
by fitting the (non-analytic part of the) Regge free energy with the optimal 
parameters in Eq. pTt with the above fmiction (|27p from Pi^ft — P* — A/3 to 

Pr^ght = P* + ar^Q 

1.797 < ci < 1.968 , C2 = -0.185 , 0.106 <a< 0.122 (28) 

the agreement appears to be very good when compared with the recent estimate 
in [55] in which the authors found Uobs ~ 0.114(6) so that 

K,.-a| ^^ (29) 
aobs 100 ^ ^ 



Figure 3: Bold points represent a sampling of Kallen-Lehman free energy 
^31) \P) while the continuous line corresponds to the critical part of the 
expected free energy Fcrit versus /3. 




Further theoretical as well as experimental determinations of a can be found 
in [P : the deviation Aa of our prediction (|28p from observations appears to be 
less than 7% in all the more recent values. The following figure Fig. ([3]) in which 
we draw both the graph of the (non-analytic part of the) Regge free energy 
and of Fcrit confirm that in the parts of the graph of Fcrit in which the 
dependence on a is important the agreement is satisfactory. 

^The fit result around the critical point is stable when 0.005 < A/3 < 0.006. In this range 
the agreement doesn't get worse than in eq. II29I I . 
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6 Further possible improvements 



The first obvious improvement is to keep all the Regge exponents C3, C2; Ci? ^-iid. 

without fixing a priori any of them (while here, for the sake of simplicity, 
we set both and C,^ to 1/2, the 2D Ising values). As it has been already 
stressed, it is also natural to explore different patterns of duality breakings. 
One can keep in the parameters space the parameter do without fixing it to 1 
(the 2D Ising value) . In the ansatz ^ the terms corresponding to the 2D Ising 
model ([2]) (which, in the Kallen-Lehman view comes from the single particles 
discrete part of the spectrum [5]) has a natural interpretation as a "boundary 
terms" . In other words, the 2D-like term with a single integral has its origin 
in the boundary contributions in the same way as in the 2D case, the trivial 
one-dimensional term has its origin in the boundary. It is conceivable that the 
patterns of duality breakings in the bulk and at the boundary could be different. 
This would lead to a Kallen-Lehman free energy in which the function /cg/ / [P) 
(see Eqs. ([TT|). (|15l) and (HH])) appearing in the purely three dimensional term 
(the double integral) has a c?o different from the k^f f (/?) appearing in the two 
dimensional term (the single integral in Eq. Two different keff in the 

single and in the double integrals (each one with its own do) could describe 
a situation in which the pattern of duality breakings in the bulk and in the 
boundary are differeniH. We expect that more general duality breakings could 
improve of two or three order of magnitude the (already quite good) agreement 
at low temperatures. As a matter of fact, all these natural ways to improve 
the numerical results (which add at least two parameters to the computations) 
would need a big cluster of CPUs since without a cluster the CPU's time needed 
would be too long. 



7 Conclusions and perspectives 

In this paper a Kallen-Lehman approach to 3D Ising model has been investigated 
numerically in the realm of a minimal duality breaking. It has been shown that 
one can fix three parameters of the model to get an excellent agreement at 
high temperatures. With the same set of parameters, the agreement at low 
temperatures appears to be very good. Remarkably enough, with the same set 
of parameters, the predicted critical exponent a has a relative deviation with 
respect to the most recent determinations of the order of the 7%. We believe 
that the present results provide one with a strong evidence that the application 
of the present methods to the study of the 3D Ising model is very promising. To 
the best of our knowledge, there are no other analytical methods able to give 
reliable informations at the same time at high temperatures, at low temperatures 
and at the critical point. More general patterns of duality breakings are worth 
to be investigated since they would further improve the already very satisfactory 
agreement with MonteCarlo data. 

*Also worth to be investigated is the case in which do has a smooth dependence on the 
temperature in order to describe a pattern of duality breaking which changes with /3. 
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